Limit of Spin Squeezing in Finite Temperature Bose-Einstein Condensates 



O 

(N 

< 



C/2 

a 

o 



oo 



o 



X 



A. Sinatra/ E. Witkowska,^ J.-C. Dornstetter/ Li Yun/ and Y. Castin^ 

' Laboratotre Kastler Brossel, Ecole Normale Superieure, UPMC and CNRS, Pans, France 
^Institute of Physics, Polish Academy of Sciences, Warszawa, Poland 

We show that, at finite temperature, the maximum spin squeezing achievable using interactions 
in Bose-Einstein condensates has a finite limit when the atom number N ^ oo ai fixed density and 
interaction strength. We calculate the limit of the squeezing parameter for a spatially homogeneous 
system and show that it is bounded from above by the initial non-condensed fraction. 
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Atomic clocks based on cold alkali atoms in two hy- 
perfine states a and b are widely used as frequency stan- 
dards. When atoms in uncorrelated quantum states are 
used, the clock precision is limited by the so-called pro- 
jection noise, resulting from the quantum nature of the 
collective spin S, i.e. the sum of the effective spin 1/2 of 
each atom. This limit is actually already reached in most 
precise clocks [1] . Spin squeezing 2] amounts to creating 
quantum correlations among the atoms so as to increase 
the precision of the atomic clock beyond this standard 
quantum limit. The relative improvement on the vari- 
ance of the measured frequency Aw^j, defines the spin 
squeezing parameter Q . Spin squeezing in atomic en- 
sembles was first obtained by quantum non-demolition 
measurements 0,11. Recently a significant amount of 
spin squeezing (e.g. 6 or 8 dB) has been achieved using 
atoms in a resonant optical cavity Q or exploiting atomic 
interactions in bimodal Bose-Einstein condensates 
The ultimate limits of the different paths to spin squeez- 
ing are still an open question. We determine here the 
influence of the non-condensed fraction for spin squeez- 
ing schemes using Bose-Einstein condensates |8l4lO|. 

A central issue is the scaling of the squeezing for large 
atom numbers. Most studies are based on a two-mode 
description [2] . In this case the squeezing parameter op- 
timized over time Cbest tends to zero (infinite metrology 
gain) for TV — oo as S,^^^^ ~ iV~^/^. The first analysis 
of squeezing at finite temperature fllj used a large N 
expansion in a Bogoliubov-like approach and could not 
predict any deviation of spin squeezing from the two- 
mode model. Here, using fully non-perturbative semi- 
classical field simulations and a powerful formulation of 
Bogoliubov theory in terms of the time dependent con- 
densate phase operator [l^ . we find on the contrary a 
dramatic effect of the multimode nature of the field: For 
a spatially homogeneous system in the thermodynamic 
limit, the two-mode scaling ^^j,^^ ^ JSf-'^/^ turns out to 
be completely irrelevant, and the spin squeezing has a 
finite optimal value that we determine analytically. 

The physical problem that we face is the dynamical 
evolution of a finite temperature Bose condensed gas af- 
ter a pulse tt/2 that puts each atom in a coherent super- 
position of two internal states a and b. This produces a 
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FIG. 1: (Color online) Spin squeezing parameter as a func- 
tion of time after the pulse mixing the states a and b, for 
N — 10^ Rb atoms (s-wave scattering length a — 5.3 nm), 
in a harmonic trap with oscillation frequency uj/2tv = SOHz. 
The Thomas-Fermi chemical potential is fi — 15.36/kiJ. Finite 
temperature semi-classical field simulations with initial non- 
condensed fractions: (a) {Nnc)/N — 0.34 (black solid line), 
(b) {Nnc)/N = 0.20 (red line), (c) {Nnc)/N = 0.09 (blue 
line), corresponding to ksT/^ — 2.08, 1.53, 1.17 respectively. 
Dashed line: Two-mode theory for comparison. 



non-equilibrium state that has a non-trivial evolution due 
to the atomic interactions inside each internal state. For 
simplicity, it is assumed that there is no cross-interaction 
between a and b atoms. 

Semi- classical field simulations - In Fig[l]we compare 
the two-mode theory with semi-classical field simulations 
at finite temperature in a trap. The gas is initially in 
state a at thermal equilibrium. In that state, we assume 
that thermal fluctuations dominate over quantum fluc- 
tuations and we use a classical field description isl- 16j 
with an energy cut-off at kgT. The initial field tpa 



then 



randomly samples the thermal equilibrium classical field 
distribution for the canonical ensemble at temperature 
T. For the initially empty state b, inspired by the trun- 
cated Wigner approach 17, we represent the vacuum 
by a classical field V'h^' having in each mode independent 
Gaussian complex fluctuations of zero mean and variance 
1/2. A sudden tt/2 pulse mixes the initial fields ipa''' ^nd 



■i/'b^'* so that, at time t = 0^ 



i.e. just after the pulse: 



V'a,.(0+) = -^[<lTA.a 



(0)i 



(1) 
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At later times, each field evolves independently according 
to the non-linear Schrodinger equation 
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(2) 

This corresponds to a harmonically trapped gas with 
same oscillation frequency u and same coupling con- 
stant g — ATrh?a/m for the two internal states, where 
a is the s-wave scattering length. As shown in Figll] 
the squeezing is created dynamically by the interactions. 
However, even for a moderate non-condensed fraction 
(iV„c)/iV = 0.09, the best in the muhimode theory 
is larger by more than one order of magnitude than in 
the two-mode theory. 

In order to isolate the effect of the non-condensed frac- 
tion from other dynamical effects taking place in the 
trapped system, as for example the spatial dynamics of 
the condensate wave function , and to develop an 

analytical theory, we consider from now-on the homoge- 
neous case. We first use the semi-classical field model 
that has the advantage that it can be simulated exactly, 
and we generalize the results to the case of a quantum 
field in the end. The real space is discretized on a lattice 
with unit cell of volume dV, within a volume V with pe- 
riodic boundary conditions 12']. The Hamiltonian after 
the pulse for component a (and similarly for b) reads: 



(3) 



The fields have Poisson brackets ih{tpf^{r),ilj*{r')} — 
Srr' 5 fiv / dV with ^,1^ = a or b, and ak(&k) is the am- 
plitude of ijja{b) the over the plane wave of momentum k. 
In terms of the fields, the collective spin components are 



Sx+iSy = / d3r^*(r)^,(r), 



(4) 



5, = -y d'r[r,ir)Mr)-rtir)Mr)]- (5) 

The spin squeezing parameter is equivalent to the min- 
imal variance of the spin orthogonally to its mean direc- 
tion, divided by the mean spin length squared and suit- 
ably normalized. Here the mean spin is along x so that 



) + {S!)^\{{Sy+^S.^' 



(6) 
(7) 



As a first step, we performed semi-classical field simu- 
lations for different temperatures and increasing system 
sizes ,21;]. The result (not shown) is that ^^^^^ converges 
to a finite value at the thermodynamic limit: N^oo, 
V—>-oo, p, T=constant, where p = N/V is the total 
density. Five independent physical parameters are in 



the model, h/m, g or a, fcsT, N and V . From dimensional 
analysis, ^^^g^. is a function of the three independent di- 



mensionless quantities that one can form, 



and 



ksT/ pg. The existence of a thermodynamic limit then 
implies 



best 



(8) 



As a second step, we performed simulations increasing 
the density in the weakly interacting limit [l^], p — oo, 
(7 — >■ with r, pg^constant. We find that, for a given 

>bcst "^"^^^ ov^aico ao X / p ^ ^ pa^ 



ksT/pg, then scales as l/pcc \/pa^- This implies: 



best/ 



F{kBT/pg) 



(9) 



In Fig [5] we show the universal behavior ©. The cir- 
cles and the squares correspond to two different values of 
\fpc? in simulations. 

Semi-classical field analytics - We now develop an ana- 
lytical theory to explain these results. We split the fields 
after the pulse as ipa = + i^a± and similarly for 4'b- 
We introduce the modulus and phase conjugate variables 
for the condensate modes 



ao 



'Nao , bo = e" 



(10) 



and we introduce number conserving non-condensed 
fields Aa and Af, that we expand over Bogoliubov 
modes with amplitudes Cak and c^k respectively (23j : 



A„ 



'al. = (^fc^'-k + ^feC:_k) -J= (11) 



Uk + Vk 



Ek 



1/4 



2 1,2 



(12) 



Ek + pg J 2m 
The spin raising component 5*+ ~ Sx + iSy is given by 



VNaoNbo + J d^AlAi 



(13) 



Our strategy is to perform a double expansion of (5+). 
We will need terms up to ~ A'^ in the thermodynamic 
limit and up to order one in the non-condensed fraction 
(Anc)/A^- In this framework, we can approximate {Sx{t)) 
in the denominator of ([7]) by its value at t = 0+ so that 



_L,min ' 



(14) 



We sketch the main steps. In the Bogoliubov limit, the 
condensate phases at t > obey [l3| 

Oa-Ob^ {0a - eb){Q+) - [{Na ~ N^) + S\ (15) 



(0) 



(^a-eb)(0+) = -2^^^+O(A^-i) (16) 
S = J2 + (l^-kP - M) ■ (17) 
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FIG. 2: Best squeezing ^bost divided by 



as a function 



of ksT/pg. Symbols: semi-classical field simulations with 
(filled squares) and 1.94 x 10"^ (disks). 
The thermodynamic limit is reached already for = 3 x 
10** except for the lowest value of ksT/pg. Lower solid line: 



analytical semi-classical field result (|21|l . Upper solid line: 
quantum result (|25p . Inset: quantum ^bcst (solid line) and 
non-condensed fraction {Nnc)/N (dashed line), both divided 
as functions of ksT/pg. 



S is the multimode part of the relative phase derivative 
that is absent in the two-mode theory. In thermodynamic 
limit Oa ~ 8b ^ ^/VN and it is sufScient to expand the 
exponential in (|13p to second order. In the modulus of 
5+ we expand: 



VNaoNbO 



1 

2 

(0)|2 



d^r {\Aa\ 



I Ah 



(18) 



with A^tot = N + ^jj. I is the total atom number in 
the semi-classical field picture . 

Best squeezing - For the calculation of the best squeez- 
ing, one looks at the asymptotic behavior of for 
t ^ oo. One finds: 



FIG. 3: "Close to best" squeezing time tr, for i] = 0.2. Filled 
squares and disks: simulations as in Fig[21 Line: analytical 
prediction (|22p . Squares and circles: thermalization times in 
the simulations extracted from the decay of the contrast {Sx)- 



analytical approach. However, the numerical squeezing 
curve as a function of time is indeed quite flat around 
its minimum, so that it suffices in practice to deter- 
mine the "close to best" squeezing time defined as 
i^{tn) = (1 + ??)$bcst' where ry > 0. Then, according to 



T9)l . t^i is finite and given by 



2 

best 



(22) 



The "close to best" squeezing time tn ((22|) for -q = 0.2 is 
shown in Figl3]and compared to simulations. 

A last important issue is that of thermalization, ne- 
glected in Bogoliubov theory and in our analytical treat- 
ment, but fully included in the semi-classical field simu- 
lations. Indeed it is possible to reach = (1 + ri)(^^^^ 
with ^^ggj given by (|2T|) only if t^ given by (|22|) is shorter 
than the thermalization time 



with the best squeezing 



1 + 



(19) 



(20) 



which remarkably only involves the multimode part p7)) 
of the phase difference. An explicit calculation gives 
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2p 
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(solid line in FigH]). In dH]), Sk = Uk + Vk given in ([T^. 
and s^i^'^ is the equivalent quantity before the pulse ob- 
tained by replacing pg with 2pg in ([T2|) : nf^ — ksT/ef'' 
are the equilibrium occupation numbers of Bogoliubov 
modes before the pulse with e^°^ — [Ek{Ek + 2pg)]^/'^ . 

Squeezing time - From (|19p . the best squeezing is 
reached in an infinite time, which is a limitation of the 



tjj <C tthc 



(23) 



In FigU] we show the squeezing parameter and con- 
trast (Sx) across the thermalization process that brings 
the system back to equilibrium after the pulse. For the 
squeezing, we compare the simulation with (i) the full Bo- 
goliubov theory (without the analytic expansions) that 
we implement numerically for a finite size system and 
(ii) a Bogoliubov ergodic model jl2] where the ampli- 
tudes Cak, Cf,k in (jlSp sample microcanonical distribu- 
tions with number of particles and the energy {Na,Ea} 
{{Nb; Eb}) fixed to a random value set by the pulse. Note 
that the simulation agrees with the Bogoliubov model at 
short times (included the "close to best" squeezing time) 
and then converges towards the ergodic model. We ex- 
tract a thermalization time ithcrm from the contrast. As 
thermalization occurs the excited modes dephase and 



(24) 
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iV — > cx) at fixed density and interaction strength, 
contrarily to the vanishing prediction of the two-mode 
model. We find that S,^^^^ is the product of \/~pa^ and 
of a universal function of koT I pg that we calculated 
analytically, and is bounded from above by the initial 
non-condensed fraction. Our analytical treatment is re- 
stricted to evolution times smaller than the thermaliza- 
tion time, but this is enough to access ^^^.^^ as we showed 
by semi-classical field simulations (that include thermal- 
ization) over a wide range of parameters. 
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FIG. 4: (Color online) Top: Spin squeezing ^ as a func- 
tion of time. Black solid line: simulation. Blue dashed line: 
Bogoliubov theory. Red dash-dotted line: ergodic model. 
Bottom: contrast {S^) (solid line) and condensate contrast 
Re(6oao) (dashed line) as functions of time. A'^ = 3 x 10^, 



ksT/pg = 3.16, 



1.32 X 10" 



The longer time scale for the decay of (bQao) is set by 
phase spreading due to partition noise (24] plus thermal 
corrections. 

Quantum field - All our analytic calculations for the 
semi-classical field can be generalized to the quantum 
field. In particular, (I^H)) and ([2^ are unchanged and 
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At zero temperature we get 
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^2 T= 
"Sbest 



-y2--ln(V2 
6 2 ^ 



1) - TT 



0.02344 . 



In practice pa^ < 10 ^ in present squeezing experiments 
so that (HH) predicts Cbe^f" ~ 2.10"^. This value is very 
low, in particular below the limit given by particle losses 
[25} . Asymptotically for fc^T pg, ^^gj,^ identifies with 
the initial non-condensed fraction. An interesting result 
is that at any temperature the initial non-condensed frac- 
tion is larger than ^^esf these two quantities in the in- 
set of FigH Aheady for keT/pg = 2, ^^^^^ and {N^c)/N 
are within a factor three. A similar conclusion seems to 
hold in a trap, see FigH] 

In conclusion we have shown that a realistic descrip- 
tion of the limits of spin squeezing in interacting Bose- 
Einstein condensates has to be multimode: The best 
achievable spin squeezing ^^^st admits a finite limit for 
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